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^ (57) Abstract: Disclosed are apparatus and methods for delivery of transcranial magnetic stimulation. The apparatus includes a 
2? TMS coil which when energized generates an electric field substantially parallel to a long axis of the coil and substantially normal 
to a surface of the coil. Furthermore disclosed an apparatus for delivery of TMS in which a coil is adapted to a robotic member 
Q for computer-aided control and delivery. Further disclosed are methods of TMS planning and delivery in which subject images are 
£^ utilized to plan, position and orient the TMS coil for precise delivery. Disclosed also are TMS coils having unique designs to better 
^ focus and direct magnetic stimulation. 
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APPARATUS AND METHODS FOR DELIVERY OF TRANSCRANIAL 
MAGNETIC STIMULATION 

BACKGROUND OF THE INVENTION 
Field of the Invention 

5 The present invention is directed to apparatus and methods particularly suitable 

for precise aiming and delivery of magnetic stimulation, and more specifically, 
transcranial magnetic stimulation. 
Description of the Related Art 

Transcranial magnetic stimulation ("TMS") is a means of repetitively stimulating 

10 the human brain through an intact scalp and skull, i.e., non-invasively. TMS is delivered 
by passing a brief (200 microsecond), strong (10,000 volts, 6,000 amps) electrical current 
through a coil of wire (a TMS stimulator) placed adjacent to the head. The passage of 
electrical current induces a strong (2 Tesla) magnetic field which, in turn, induces 
electrical currents in nearby tissues. In the case of nerve cells, if the induced current is 

15 sufficiently intense and properly oriented, it will result in synchronized depolarization of 
a localized group of neurons (i.e., neuronal "firing")- Initially, magnetic stimulation was 
used only for peripheral nerves, in which instance it is affecting nerve fibers rather than 
neuronal cell bodies. More recently (Barker et al., 1985), magnetic stimulation has shown 
to be able to depolarize neurons in the brain. The cellular element of the brain being 

20 affected by TMS was assumed, but not proven, to be fibers rather than neuronal cell 
bodies. 

TMS has several present and potential applications, in the domains of basic 
neuroscience research and of the treatment of brain disorders. Applications for 
neuroscience research include, for example: imaging brain connectivity (e.g., Fox et al., 

25 1997); establishing inter-regional and inter-hemispheric conduction times (e.g., Meyer et 
al., 1995); testing the function of specific brain areas by means of transient functional 
disruptions, so-called "virtual lesions" (e.g., Shipley & Zeki, 1995); and, studying the 
modification of synaptic efficacy induced by repetitive stimulation, termed LTP (long- 
term potentiation) and LTD (long-term depression). Potential clinical applications 

30 include, for example: pre-operative mapping, e.g., of language related brain areas 
(Epstein (et al., 1996); testing for neuronal conduction delays due to dysmyelinating 
disorders; and, treating brain disorders by selective modification (up or down regulation) 
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of the synaptic efficacy of pathways (i.e., by inducing LTP and LTD; Wang, Wang and 
Scheich, 1996). 

At present, TMS delivery is crude. The TMS effector or stimulator (commonly 
referred to as a "TMS coil")is a wire-wound coil, typically shaped like a "B." The B- 
5 shaped coil is placed against the scalp and held in place by a human operator. For the 
primary motor cortex and primary visual cortex (small sections of the total brain surface), 
proper positioning is established by the elicited response: muscle contractions when 
stimulating the primary motor cortex; illusory lights (phosphenes) when stimulating the 
primary visual cortex. In both of these areas, the effects are very sensitive to coil position 

10 and orientation. 

For brain regions in which proper positioning cannot be determined by the 
induced effects (i.e., muscle contractions or subjective experience), position is generally 
determined by reference to a traditional pattern used for placement of EEG electrodes 
(10/20 system). The 10/20 system is based on scalp/skull landmarks which do not bear a 

15 reliable relationship to the functional anatomy of the brain. Further, when using the 10/20 
system, there has been no strategy enunciated for determining proper orientation of the 
coil. Thus, a reliable method for determining the proper position and orientation of TMS 
coil placement for brain areas lacking immediately observable feedback is needed. 
Application of TMS during radionuclide imaging (using positron-emission 

20 tomography ("PET') or single photon emission tomography ("SPECT")) has two 

important uses. First, radionuclide imaging can be used to monitor the induced response, 
determining precise location and quantifying response magnitude. This is extremely 
important for testing aiming algorithms and for determining the effect of stimulation 
parameters, such as intensity, rate, duration and the like. Second, an important use of 

25 TMS is to map brain connectivity using radionuclide imaging. For both of these 

applications, hand-held TMS delivery is inappropriate, for at least three reasons. First, 
hand-held delivery is unsafe, unnecessarily exposing the experimentor to the radiation 
used for imaging. Second, hand-held delivery is positionally unstable, degrading image 
quality by small movements of the holder. Third, hand-held delivery is intrinsically 

30 inaccurate and imprecise. 

Further, current coil designs for delivery of TMS have been mainly intuitive and 
somewhat crude. Typical coil designs consist of two loop figure eight type coils, for 
peripheral nerve and brain stimulation, four loop coils for peripheral nerve stimulation, 
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and variations in angles of these. While attention is paid to coil inductance, it is only for 
simple circuits that this may be easily calculated. 

The target field method has been used to produce minimum inductance cylindrical 
gradient coils for MRI (Turner, 1986) and has been adapted for bi-planar coils (Martens 
5 et al., 1991). Minimum power designs have also been presented (Bowtell et al.). 
However, such design methods have not been applied to the design of magnetic 
stimulation coils. 

Various combinations of circular or rectangular coil shapes have been designed. 
Figure eight type designs appear to be the most common. Further, B-shaped and slinky- 

10 type coils have also been designed (Cadwell; Lin et al.). Sections of toroids (Davey, 

Epstein, Carbunaru) with magnetically permeable cores are also known, and appear to be 
an efficient design. However, none of these provide an extremely focused field 
penetration. Two and four wing coils with a straight section joined with curves for 
peripheral nerve stimulation have been designed. (Ruohonen et al.) However, these 

15 designs are largely limited by intuition. 

Thus, fundamental limitations on the utility of TMS for research and treatment 
include a lack of methods for precise, automated aiming (positioning and orienting) and 
safe, rigid (i.e., non-human) holding of the TMS stimulator, as well as the poor suitability 
of present coils for TMS. 

20 SUMMARY OF THE INVENTION 

The present invention enhances the precision and ease with which TMS may be 
used for the diagnosis and treatment of neurological and psychiatric disorders and for 
neuroscience research. In certain embodiments, these benefits may be accomplished via 
use of specifically shaped TMS stimulators having certain properties. In certain 

25 embodiments, a robot, such as a neurosurgical robot, may be used to deliver TMS. The 
present invention also includes algorithms for treatment planning and treatment delivery, 
including: algorithms for rapidly modeling the 3-D electric field created in the brain by a 
TMS coil at any external location; cortical surface modeling (extraction and 
visualization); scalar product (electric-field vector times cortical-surface vector) 

30 computation and visualization; and merging of functional images, structural images and 
treatment-planning models (surfaces & fields). 

The present invention further includes treatment-delivery tools such as frameless 
registration of head, brain image, and robot; fully automated robotic positioning of the 
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TMS coil; robotic sensing of TMS orientation (about a manually operated tool-rotation 
axis); hardware extensions including a passive digitizing arm; a TMS tool mount; a 
passive tool-rotation axis with an orientation sensor; and a general-purpose mobHe cart. 

Additionally, inventive methods enhance the precision and ease with which TMS 
5 may be used for neuroscience research and for the diagnosis and treatment of 
neurological and psychiatric disorders. 

The inventors have determined thai the biological efficacy of transcranial 
magnetic stimulation applied to cerebrum can be estimated at any point as the scalar 
product of the induced electrical field (E, a vector) and a unit vector aligned parallel to 
1 0 the cortical columns. The unit vector is estimated as a normal (i.e., perpendicular) the 
cortical surface, as this is the known orientation of cortical columns. The biological 
efficacy of a TMS E field, then, is calculated using equation 1. 

Biological efficacy = |E| cos 0 [1 ] 

where G is the angle between E and the unit normal vector. This Cortical Column 
1 5 Cosine Aiming Principle ("CCCAP" or "CAPS" or "aiming principle") is based on the 
inventors' determination that the cortical column is the biological unit of the brain with 
the lowest threshold for TMS excitation and the well-established neuro-anatomical 
principle that the cortical columns are oriented at a right angle to the cortical surface. 
Thus, maximum biological efficacy of a cortical region of interest occurs where the 
20 induced E field is parallel to the direction of the cortical columns, i.e., normal to the 
cortical surface. 

The CCCAP includes the following principles: (1) surface grey matter (the 
cortex) is preferentially or exclusively activated; axons in the sub-cortical white matter 
are minimally activated by the TMS-induced E- field, but will conduct action potentials 

25 initiated in cortex by TMS; (2) cortical grey matter will be most effectively activated by 
an E-field oriented parallel to the columnar organization of the cortex; (3) the response 
magnitude at any cortical location is a function of the magnitude of the E-field parallel to 
the cortical columns; (4) cortex is preferentially (but not exclusively) activated by 
orthodromic E-fields, passing from the pial surface, through the soma, to the sub-cortex; 

30 antidromic E-fields, passing toward the pial surface will be less effective but not 

ineffective; and (5) the orientation of the cortical columns is macroscopically estimated as 
the normal to the true cortical surface. Of the five aiming principles, the first two 
(particularly the second) are strongly at odds with current opinion and practice in the 
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TMS community. It is to be understood that a true cortical surface is one derived from an 
anatomical image with sufficient spatial resolution and image contrast to define the cortical- 
subcortical or cortical-CSF border. Simplified (e.g., spherical) models or generalized models 
of the cortex which do not define the true cortical surface may not accurately model the 
5 orientation of the cortical columns. 

FIG. 1 is a cross-sectional view of a subject's head 5 with a B-shaped coil 10 
positioned thereabove. As shown, the head includes the scalp and skull 5, the cerebral cortex 
or grey matter 15, cerebrospinal fluid (CSF) 20, and white matter 25. Where the cortex is 
concave, folding inwardly away from the scalp and skull, it is termed a sulcus 35 (pi. sucli). 

10 For present modelng purposes, the most important components of the cortex are the vertical 
neurons 40, which are oriented perpendicular to the brain's cortical surface (i.e., 
perpendicular to the interface between cortex 15 and CSF 20). The vertical neurons 
collectively form the cortical columns (not shown), which are the dominant anatomical and 
physiological features of cerebral cortex (at the microscopic level), being present in all 

1 5 regions of cortex in all mamallian species. (Only a few vertical neurons are illustrated in 
Figure 1). 

In FIG. 1, the TMS coil 10 is positioned and oriented to create an induced E-field that 
is perpendicular to the brain cortical surface (and parallel to the vertical neurons) in the 
sulcus (C), but parallel to the cortical surface (and perpendicular to the vertical neurons) at 

20 the crown of the gyrus (A). The magnitude of the E field is weaker at C than at A, because 
the distance from the surface of the TMS coil 10 is greater at C than at A. Relative to the 
orientation of the vertical neurons forming the cortical columns the E field can be 
decomposed into vertical (Ey) and horizontal (Eh) components, which are parallel to the 
vertical neurons and horizontal fibers, respectively. By Principal (2), activation at any 

25 cortical site is a function of Ey with negligible E h effect. Thus, E has no effect (no Ey 

component) where E is perpendicular to the column (A); intermediate effects for intermediate 
relative orientations (B); and maximal effect where E is parallel to the column (C). 

Prior art (Brasil-Neto et. al 1992; Mills et al., 1992) has empirically demonstrated that 
coil orientation computed in the manner just described is optimal for one brain area (the 

30 primary motor cortex), but provided a rationale not applicable to brain regions other than 
primary motor cortex and inconsistent with the above explanation. 
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In an example embodiment, the CCCAP may be used to allow the cortical 
excitation effects of TMS to be computed in advance for any position and orientation, 
thereby enabling computer-aided aiming of TMS. However, in other embodiments, the 
CCCAP may be used for manual aiming and orienting. The CCCAP may also be used to 
5 normalize (correct) observed biological effects for the angle of intersection with the 
cortex, when aiming is done conventionally (i.e., not with CCCAP) but images are 
obtained showing the relationship between the TMS coil and the cortical surface. 

In an example method, image-guided, computer-aided implementation of the 
CCCAP for TMS delivery may be effected in accordance with the following steps. First, 

10 an imaging stage is performed, in which a high-resolution, anatomical image (e.g., a 3D 
Tl-weighted image) of the subject's head is obtained. Preferably, sufficient grey-white 
contrast is obtained to allow detailed cortical-surface segmentation. Next, a functional 
image (either functional MR] or PET, for example) data set is obtained during conditions 
(task/control pair) which selectively activate the cortical region of interest (e.g., repetitive 

1 5 hand movement to activate the supplementary motor area). 

Next, a modeling stage is performed to identify surfaces within the anatomical 
and functional images, to be used for registration of the functional and anatomical images 
to one another and, subsequently, to the patient. In the anatomical image (e.g., an 
anatomical MRI), the scalp surface and the brain's cortical surface are segmented and 

20 modeled as polygon-mesh surfaces. The brain-surface is modeled at high-resolution, to 
provide an accurate, detailed representation of the interface between cerebrospinal fluid 
(CSF) and cortical grey matter, as this is critical for establishing the orientation of the 
cortical columns. Anatomical surface extraction and modeling is done in a manner 
keeping both surface models (scalp and brain) in register with the original image and, 

25 thereby, with one another. In the functional image (either PET or fMRI), the brain 
surface is segmented and modeled as a polygon mesh surface (but with less detail than 
the model derived from the anatomical image) and the targeted site is identified. The two 
brain surface models (anatomical-image-derived and functional-image-derived) are co- 
registered, thereby bringing the target site in register with the two MRI-derived cortical 

30 surface models: brain and scalp. This comprises a conjoined functional/anatomical 
model. 

Then, a 3-D model of the TMS-coil's physical surface and the 3-D E-field induced 
by the TMS coil is created, and this TMS coil-surface/E-field model is superimposed on 
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the conjoined functional/anatomical model. The TMS coil-surface/E- field model is 
positioned and oriented so as to obtain maximum biological efficacy (as defined by the 
CCCAP) at the target point, while keeping the coil surface outside the scalp surface. The 
position and orientation of the TMS coil-surface relative to the conjoined 
5 functional/anatomical model may then be stored for subsequent use. 

In an example embodiment, this data may be used to perform TMS delivery in 
accordance with the following. First, the subject is placed in the treatment position and 
the head immobilized, for example, using a thermoplastic mask (Fox et al., 1985). With 
the subject in the treatment position, a 3-D digitizer (e.g., a passive arm digitizer) may be 

10 used to collect a series of points on the scalp surface. These points are used to create a 
model of the scalp surface in the subject's current head position. This model is registered 
to the conjoined functional/anatomical model previously created, using a rapid, surface 
matching algorithm, such as a convex hull algorithm (Lancaster et al., 1999). The manual 
digitizer may again be used to collect a set of specific reference points on the surface of 

15 the TMS coil, which may be mounted on a multi-joint, calibrated armature, either passive 
or robotic. 

As the optimal position for the TMS to achieve supra-threshold stimulation of the 
target location has been previously computed, the translations and rotations needed to 
move from the present position to the desired position are computed. If the TMS is held 
20 by a passive armature, movement is executed manually, with the readout of the 

coordinates of the armature as feedback. If the TMS is held by a robot, movement is 
executed automatically, after a safe pathway (avoiding contact with the subject or other 
obstacles) has been computed. When properly positioned and oriented, TMS delivery is 
effected. 

25 The advantages of this method are several. First, placement is extremely precise 

(± 1 mm), much more precise than is possible using hand-held aiming. Second, placement 
can be computed in advance, rather than by trial-and-error, as is done with hand -held 
aiming for areas with a measurable behavioral response. Third, once computed, a position 
can be precisely re-established for subsequent treatments. Fourth, the entire process is 

30 mathematically speci fied and suitable for computer implementation. Fifth, the procedure 
is suitable for application to any cortical location, not just brain areas in which aiming 
accuracy can be confirmed by an overt behavioral response. Sixth, once the necessary 
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brain images are acquired, trajectories for stimulating any number of functional areas can 
be computed. 

in an example embodiment, TMS delivery may be effected by means of a robotic 
arm. Using a robotic arm to move, aim and hold the TMS provides many benefits. First, a 
5 robotic arm allows the TMS coil to be placed and oriented with accuracy and precision of 
- 1 mm (location) and — 1° (orientation), far exceeding hand-holding. This degree of 
accuracy and precision is crucial for implementing a high-precision aiming algorithm. 
Second, a robotic arm allows the TMS to be positioned automatically, i.e., under 
computer control. Third, a robotic arm allows the TMS to be positioned rapidly, saving 

1 0 time for the patient, experimentor or clinician. Fourth, a robotic arm can accurately re- 
establish its position on sequential days, which is needed for experiments and treatments 
applied over a series of sessions. Fifth, a robotic arm allows the TMS to be moved rapidly 
from one position to another, allowing treatment of more than one location in a single 
session. Sixth, a robotic arm allows the TMS to be held immobile for long periods of 

1 5 time, which a person cannot practically or comfortably achieve. Seventh, a robotic arm 
allows TMS to be delivered during PET imaging, without exposing a human holder to 
radiation. 

The present invention is also directed to an inverse design method to produce a 
desired electric field profile for transcranial/peripheral nerve magnetic stimulation. Such 

20 a method may be used to design a variety of coils. Such coils may include multiple 

winged coils having a concentrated bundle of wires at the center and smooth arcs splayed 
therefrom at an increasing distance, maintaining the minimum inductance and/or 
minimum power dissipation possible for a given field profile. Similarly, other coils may 
be used to produce a spatial gradient across the nerve, maintaining minimum inductance 

25 and/or minimum power dissipation. 

Coils may be designed in accordance with the teachings of the present invention 
to incorporate one or more of the following characteristics which aid in the delivery of 
transcranial magnetic stimulation. In certain embodiments, a coil may be designed to 
provide for minimum inductance for a given set of field constraints. Further, a coil may 

30 be designed to dissipate a minimum amount of power for a given set of field constraints. 
Similarly, a coil may be designed to provide for minimum inductance and minimum 
dissipation for a given set of field constraints. In certain embodiments, a coil may be 
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designed with negative turns to reduce the electric field on a patient. Such negative turns 
also may be provided in a separate layer to similarly reduce electric field on the patient. 
BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention may be better understood, and its numerous objects, 
5 features, and advantages made apparent to those skilled in the art by referencing the 
accompanying drawings. 

FIG. 1 is a cross sectional diagram of a subject's brain and a B-shaped TMS 
stimulator positioned thereabove having an induced electric field parallel to the short axis 
of the coil. 

10 FIG. 2 is a diagram of a typical current configuration for an O-shaped coil. 

FIG. 3 is a diagram of a typical current configuration for a B-shaped coil. 
FIG. 4 is a simplified top view of TMS coil orientations used in a prior art study 
of primary motor cortex. 

FIG. 5 is a simplified top view of TMS coil orientations used in a prior art study. 
15 FIG. 6 is a graph of the results of delivering TMS during PET applied to the foot 

area of primary motor cortex at 0° (perpendicular to cortex), 30°, 60°, and 90° (parallel to 
cortex) . 

FIG. 7 is an isometric view of an example TMS stimulator according to the 
present invention. 

20 FIG. 8 is a diagram of a TMS attached to a robot and placed behind a PET 

scanner, to acquire a PET functional imaging study during TMS stimulation. The upper 
right quadrant of the PET is shown in a cut-away view, to allow the robotic TMS system 
to be illustrated. 

FIG. 9 is an isometric view of a TMS tool holder according to the present 
25 invention. 

FIG. 1 0 is a flow diagram of an example TMS treatment planning system 
according to the present invention. 

FIG. 1 1 is a flow diagram of an example TMS treatment delivery system 
according to the present invention. 
30 FIG. 12 is a representation of a TMS coil in accordance with the present 

invention. 

FIG. 13 is a representation of a TMS coil in accordance with the present 
invention. 
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FIG. 14 is a diagram of an alternate coil embodiment according to the present 

invention. 

FIG. 1 5 is a graphical representation of the electric field magnitude/current of the 
coil of FIG. 14 along the x-axis. 
5 FIG. 16 is a graphical representation of the electric field magnitude/current of the 

coil of FIG. 14 along the y-axis. 

FIG. 17 is a diagram of an alternate coil embodiment according to the present 
invention. 

FIG. 1 8 is a close up view of the center portion of the coil shown in FIG. 1 7. 
10 The use of the same reference symbols in different drawings indicates similar or 

identical items. 

DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 

Aiming is the most fundamental, unresolved issue for scientific and medical uses 
of TMS. TMS is most often "aimed" merely by observing its behavioral effects. For 

1 5 example, primary motor cortex is frequently identified by adjusting coil 

position/orientation to achieve a contraction of the desired muscle (e.g., abductor pollicis 
brevis) at the lowest stimulus voltage. Very recently, image-guided aiming has been 
introduced (below). Even when image-guided, coil orientation relative to brain anatomy 
has either been ignored (Paus et al., 1997) or empirically derived (Krings et al, 1997), but 

20 it is not based upon a general theory of aiming. Preliminary conjectures as to the mode of 
interaction between TMS and the brain have been put forward, and emphasize 
interactions with fibers running parallel to the cortical surface (i.e., horizontal relative to 
the vertical orientation of the cortical columns). While many accept the horizontal-fiber 
hypothesis discussed (below), no thorough treatment of its implications for TMS aiming 

25 exists. In fact, its predictions do not agree with a growing number of observations. The 
status and shortcomings of TMS aiming theory and practice are synopsized here, as 
further background to the present invention. 

In a planar electromagnetic coil, the induced E-field parallels the current in the 
coil, is maximum in the plane of the coil and falls off rapidly with distance from the 

30 surface of the coil. Presently known TMS coils are made in two basic geometries: (1) 

circular (O-shaped, as shown in FIG. 2); and (2) double ring (B-shaped, as shown in FIG. 
3). In air, a circular coil induces a circular E-field. Because the current density is evenly 
distributed around the coil, the E-field is curved and maximum at the edge of the coil 
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rather than at the center. This lack of a focal linear "sweet spot" makes a circular coil 
very difficult to aim. 

A much more focal E- field can be created by placing two circular coils, current 
flowing in opposite directions (i.e., clockwise and anti-clockwise), next to each other, 
5 creating a B-shaped coil. For this coil, the E-field is enhanced where the coils are near 
each other, because the electric fields from each coil sum. The result is a focal E-field in 
the center of the coil, oriented parallel to the central region (short-axis) of the B-shaped 
coil (FIG. 3). The strong, focal E-field makes the B-shaped coil much easier to aim and 
much more likely to induce a focal brain activation (Roth et al., 1991). 

10 Early O-coil studies (Barker et al., 1 985; Barker et al, 1987) observed that motor 

responses could be obtained from either hemisphere by reversing the direction of current 
flow. They concluded that this was due preferential sensitivity of neurons to orthodromic 
currents, recognizing that neural structures "are more likely to be stimulated if they are 
oriented parallel to the electrical field lines" (Barker et al., 1987). Nevertheless, early 

1 5 researchers remained agnostic on the brain location (cortical vs. subcortical; sulcal vs. 
gyral) and specific neural elements being stimulated. For example, Yeomans (1990) 
comments that "localization and orientation effects of magnetic stimulation have not 
been explained neurally" (pg. 141). In a similar vein, Roth et al. (1991) concludes: 
"Although the use of magnetic stimulation is growing rapidly, the technique has been 

20 applied without a complete theoretical understanding of the induced electric field 
distribution"; and, "none of these studies define a definitive relationship between the 
electric field distribution in the motor cortex and the resulting transmembrane potential 
induced in the cortical neurons." 

To this day, many remain agnostic on this matter. As a consequence, TMS 

25 position (placement and orientation) follows no coherent theory. Some users position the 
TMS coil by behavioral optimization (e.g., minimum-threshold muscle twitch); some 
position carefully, but ignore orientation, using the B-nose convention (discussed below); 
some orient by reference to prior behavioral-optimization studies; none use an aiming 
theory. 

30 Day (et al. 1989) were the first to postulate a specific site/mechanism of action for 

TMS. In the context of a study of primary motor cortex, Day reasoned that "stimuli are 
likely to activate those neurons nearest the stimulating electrode, that is, those on the 
convexity of the precentral gyrus". Day knew that an O-shaped coil, positioned flat 
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against the scalp, produced a circular E-field tangent to the gyral crown and at a 90° angle 
(horizontal) to the (vertical) columnar organization of the cortex at the gyral crown. Day 
accepted the long-established principle that "a voltage gradient parallel to the long axis of 
the neuron is the most favorable" (Day et al., 1 989). Consequently, Day postulated that 
5 horizontal fibers (interneurons, pyramidal tract axon collaterals and afferent axons from 
cortical and subcortical sites) at the gyral crown must be the directly activated neural 
elements, with pyramidal neurons (projecting to the spinal cord) being secondarily 
activated. Note that Day's inferences were based on the assumption that TMS would 
only be able to activate the brain region nearest to the coil surface, i.e., the gyrus. 

10 However, the data of Fox (et al. f 1997), Paus (et al. 1997) and the EXEMPLARY 

STUDY below using PET imaging to detect TMS effects on the brain, clearly indicate 
that TMS applied with a B-shaped coil preferentially stimulates sulci, rather than gyri. 
Thus, the premise of Day's reasoning was incorrect. 

Experimental support for the horizontal-fiber hypothesis is notably lacking. Day 

15 presented no experimental evidence of horizontal fiber activation. The cortical electrical 
stimulation literature (which excites cortex by passing applying electric potentials directly 
to the scalp) has consistently utilized current applied perpendicular to the cortical surface 
at gyral crown (i.e, in the orientation prescribed by the CCCAP) and never horizontally, 
as would be predicted by Day's horizontal-fiber hypothesis. In fact, no evidence has ever 

20 been provided that horizontal-fiber activation results in pyramidal neuron activation or, 
more importantly for the present proposal, in columnar excitation comparable to that of 
physiological activation or vertical electrical cortical stimulation. 

Theoretical flaws in the horizontal-fiber hypothesis are readily identified. The 
horizontal-fiber hypothesis postulates that pyramidal motor neurons, located on the 

25 anterior bank of the central sulcus, are activated by means of horizontal fibers on the 

crown of the pre-central gyrus. This is a distance of 0.5-2.0 centimeters. Yet, horizontal 
fibers extend only 1-2 mm (Jones and Wise, 1978; Jones, 1981). Further, horizontal 
fibers are isotropic, extending uniformly in all directions within a plane parallel to the 
cortical surface. The isotropism (lack of a preferred orientation) of the horizontal fibers 

30 should translate into a lack of a preferred orientation for TMS, as the E-field should 

excite a roughly equal fraction of the total horizontal-fiber population in any orientation. 
This has been clearly disproven by some of the very studies which espouse the 
horizontal-fiber hypothesis. Finally, the horizontal fiber hypothesis is weakened by the 
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fact that the predominant horizontal element providing for horizontal interactions is the 
basket cell, which "should be considered a class of large inhibitory interneuron" (Jones, 
1981). Thus, the hypothesis that TMS excites pyramidal cells by means of horizontal 
fibers is quite unlikely. 
5 Current common practice is to place a B-shaped coil flat against the scalp 

approximately over the area to be stimulated, as shown in FIG. 1 . Coil orientation is 
stereotyped, taking no account of the potential interactions between cortical geometry and 
E-field orientation. Specifically, the TMS coil is customarily positioned such that the 
short-axis of the coil (flat side of the B) faces the nose and the handle of the coil (round 

10 side of the B) faces the occiput: the so-called "B-nose" position. The logic of this 
orientation convention is not stated in the literature. It may well originate from the 
physical convenience of standing behind the subject with the coil handle facing the 
experimenter and the cable draped over the experimenter's shoulder. 

Hand-held, non-image-guided application of TMS is the norm. This rather casual 

15 approach to stimulation position and orientation likely has several contributing factors. 
First, hand-held aiming is simple to implement and inexpensive. Second, many users of 
TMS conceive its effects to be rather diffuse. Consequently, they see little need to be 
highly precise in aiming. Direct evidence that TMS' effects are highly localized only 
appeared recently, by PET imaging during TMS (Paus et al., 1 997; Fox et al., 1997). 

20 Third, no general aiming principle had been enunciated to drive greater precision. 
Fourth, no systems for facilitating precise, image-guided aiming are available. 

Paus (et al. 1997) used a modified commercially available passive digitizing arm, 
namely a VIEWING WAND® (ISG Technologies, Toronto, Canada) to perform a first- 
generation form of image-guided TMS. TMS was delivered to the average location of the 

25 frontal eye fields (FEF), defined stereotactically (i.e., in standardized coordinates) as the 
mean of eight, previously published group-mean, PET-activation studies (Paus et al., 
1996). Cortical surface geometry (sulcal/gyral location and orientation) was not taken 
into consideration either for placement or for orientation. Orientation was "B-nose". 
Anatomical MRI was obtained but was used solely to establish the inverse transformation 

30 from the standardized space (Talairach and Tournoux, 1988) to the individual's native 
image space. No functional markers (e.g., PET during eye movements) were obtained in 
the subjects either prior to or during the TMS study. The response location was imaged 
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wilh PET, but was not interpreted relative to group or individual cortical geometry nor 

was any discussion of an aiming theory included. 

Krings (et al. 1997) used a modified SURGICAL ARM® (Radionics, Burlington, 
MA) to stimulate primary motor cortex. The purpose of this study was to compare 
5 response locations for TMS and EBS, validating TMS for pre-operative mapping of 
motor cortex. Position was adjusted to achieve a consistent motor response (thumb 
twitch). Orientation was based on anatomical MRI, positioning a B-shaped coil with the 
short axis perpendicular to the central sulcus, in keeping with the recommendations of 
Brasil-Neto et al. (1992) and Mills et al. (1992). However, Krings offers no aiming 
10 theory. The orientation rules which Krings cites (Brasil-Neto et al., 1992; Mills et al., 
1992) are specific to primary motor cortex and do not constitute a generalized aiming 
theory. 

Extremely important factors for any TMS aiming system to address are alterations 
of the induced fields by the stimulation apparatus or the stimulated object. The TMS- 

1 5 induced B field falls off rapidly with the distance from the coil. The induced B field can 
be altered by nearby ferro-magnetic materials, via the creation of secondarily induced B- 
fields. The magnitude of a secondarily induced B field is a function of the \i of the metal 
and the amount of metal; the shape of the secondarily induced B-field is a function of the 
shape of the object. The secondary B field and the TMS B field will sum, which can 

20 distort the net B field in a complex manner. Materials with high \i are metals. Thus, it is 
desirable to use a robotic arm having low \x. Furthermore, interaction between coil and 
arm may be minimized by creating distance between the TMS coil and arm. Biological 
tissues, in general, are paramagnetic, having very small \i; thus, they do not significantly 
distort the TMS-induced B field. 

25 Unlike the B field, the brain's TMS-induced E field is subject to alterations from 

biological tissues, as follows. Following the TMS-induced E field, currents flow. At the 
interfaces of tissues whose conductivities differ, charge accumulates. For the head there 
are three important conductivity interfaces: 1) air and scalp; 2) scalp and skull (outer 
table); and 3) skull (inner table) and brain. While minor conductivity differences do exist 

30 among the several extracranial tissues (dermis, subdermal fat, galea) and intracranial 
tissues (i.e., meninges, blood vessels, grey matter, white matter and CSF), these 
differences are small relative to the three major interfaces above. Soft tissues (scalp and 
brain) are weak conductors; air and skull are non-conductors (insulators). Charge build up 
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occurs, therefore, at the air-scalp, scalp-skull, and skull-brain boundaries. The shape of 

the head and skull determines the geometry of the charge build up and, thereby, the 

precise geometry of the E-field distortion. The greatest change in the TMS-induced E 

field will occur nearest the accumulated charge. Research continues in this area. 

5 The CAPs predict the results of the two orientation-optimization experiments in 

the TMS literature (Brasil-Neto et al., 1992; Mills et al., 1992), and the two PET/TMS 

papers published to date (Paus et al., 1997; Fox et al., 1997), and an exemplary study 

testing the aiming principles. 

The optimal orientation for TMS excitation of primary motor cortex for hand 

10 (Ml -hand) was simultaneously reported by two independent research teams (Mills et al., 
1992; Brasil-Neto et al., 1992). Mills stimulated the left hemisphere; Brasil-Neto 
stimulated the right. Both groups used the "B-nose" orientation as the reference (0°) 
orientation, rotating the coil in 45° increments. Both groups found the optimal 
orientation to be rotated medially 45°, and interpreted this as being perpendicular to the 

15 central sulcus, as shown in FIG. 4. The orientation of the E-field (optimal current vector) 
was normal to the cortical surface of the "presumed central sulcus", as predicted by the 
CCCAP. Prior art, however, has espoused no general aiming principal nor provided an 
explanation of the observed effect consistent with known physiology or anatomy. 
Neither group interpreted this in reference to the columnar organization of cortex, with 

20 the E field being aligned with the columns of primary motor cortex, on the anterior bank 
of the central sulcus. Rather, Day's assumption of selective activation of gyral crowns 
and, therefore, the horizontal-fiber hypothesis was accepted and cited. Mills recognized 
that isotropically distributed horizontal fibers cannot explain the observed orientation 
effect and postulated that horizontal fibers might have an in-plane orientation preference 

25 at right angles to the central sulcus. In effect, he postulated a cortical "row" system 

composed of horizontal fibers, which has never been demonstrated, either anatomically or 
physiologically. 

Brasil-Neto (1 992), on the other hand, postulated that "horizontal interneurons 
which are aligned perpendicular to the central sulcus are preferentially activated by 
.30 magnetic stimuli", but provided no rationale or supporting data. Both Mills and Brasil- 
Neto cited studies which indicate the existence of horizontal fibers and that some fraction 
of these fibers are in the appropriate plane. Neither presented a rationale for preferential 
orientation (i.e., non-isotropism) of horizontal fibers nor a rationale for their selective 



-15- 



WO 02/089902 PCT/US02/14157 
orientation perpendicular to the central sulcus. Finally, neither Mills nor Brasil-Neto 
addressed the great distance ( 1 -2 cm) between the presumed site of stimulation (gyral 
crown) and the location of human motor cortex (anterior bank of central sulcus) nor the 
inhibitory nature of horizontal intemeurons. Thus, both studies are problematic for the 
5 horizontal-fiber hypothesis, even though both cited the horizontal-fiber hypothesis as 
their underlying theoretical construct. 

On the other hand, these orientation-effect results are exactly predicted by aiming 
principles 2 and 3, which predict maximal activation when the E-field is perpendicular to 
the cortex and a graded response proportionate to the total Ev. When a B-shaped coil is 

10 positioned over a sulcus, the vertical component (Ev) of E is maximal when E is 

perpendicular to the sulcus, as shown in FIG. 1. This was observed by both orientation 
studies (Mills et al, 1992; Brasil-Neto et al., 1992), as shown in FIG. 4. 

The TMS/PET experiments of Paus (et al., 1997) and Fox (et al, 1997) provide 
strong evidence for the columnar aiming principles set forth above. Paus applied TMS 

15 with a B-shaped coil to the mean location of the frontal eye fields (FEF). The resulting 
activation was recorded as an increase in cerebral blood flow using H2 15 0 PET. For 
purposes of the aiming principles, Paus* most notable finding was that the induced 
activation was sulcal (not gyral), lying approximately 3 cm deep to the gyral crown. This 
is in agreement with the prediction of principle 2 and contradicts the prediction of the 

20 horizontal fiber hypothesis. 

Fox and colleagues applied TMS with a B-shaped coil to primary motor cortex. 
The coil was positioned and oriented behaviorally, to elicit a contraction of the abductor 
pollicis brevis at the lowest possible threshold. The resulting activation was recorded as 
an increase in cerebral blood flow using H 2 !5 0 PET. Data underwent a group-mean 

25 analysis to create z-score activation images overlaid on the averaged MRI. Again, the 
induced activation (at the site of TMS application) was sulcal, lying 2-3 cm deep to the 
gyral crown. Here, the activation site was the anterior bank of the central sulcus, in 
agreement with known physiology and histology and clearly supporting principle 2. 

EXEMPLARY STUDY 
30 The inventors have carried out two studies of aiming principles 2 and 3: one 

studying E-field intensity (Experiment 1) and a second studying E-field orientation 
(Experiment 2), as follows: 
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Experiment 1 : By Principle 3, maintaining the E-field at a field orientation but 
varying the strength of the field should vary the degree of TMS-induced neuronal 
activity. Eight normal volunteers were imaged by PET while undergoing 3-Hz TMS 
stimulations across a range of intensities. Stimulation was applied to the Supplementary 
5 Motor Area on the medial surface of the right frontal lobe of the cerebral hemisphere. 
The E-field was oriented perpendicular to the cortical surface, which was also 
perpendicular to the mid-sagittal plane and exactly 90° rotated from the orientation 
required by the "B-nose" principle, described above. This stimulation produced strong 
activation on the medial-bank cortical surface at the location in which the E-field was 

1 0 most perpendicular to the cortical surface, as predicted by Principles 2 and 3. Response 
magnitude was highly correlated with TMS E-field strength, as predicted by Principle 3. 

Experiment 2: By Principles 2 and 3, varying the orientation of the TMS E-field 
relative to the cerebral cortical surface should modulate response magnitude, with the 
greatest response being observed when the E-field is perpendicular to the cortical surface. 

1 5 Twelve normal volunteers were imaged by PET while undergoing 3Hz TMS at fixed 
intensity but varying orientation: 0°, 30°, 60° and 90° relative to the cortical surface. 
Stimulation was applied to the foot area of primary motor cortex. As predicted by 
Principles 2 and 3, the response was located on the bank of the central sulcus, rather than 
on the gyral crown. Also as predicted by the same principles, response magnitude was 

20 highly correlated with TMS E-field orientation, with the greatest response occurring 
when the field was at 90° to the cortical surface, as shown in FIG. 6. 

To take advantage of the aiming principles set forth herein, it is desirable to use a 
TMS stimulator having properties which allow it to generate an E field having a 
maximum at a clearly defined location, said location being easily placed in close 

25 alignment with a subject's head, so that the E field is directed substantially parallel to the 
cortical column of interest. Furthermore, in an example embodiment, such a stimulator 
(or a standard TMS stimulator) may be used in connection with a robotic system to 
permit for accurate delivery of the E field. It is to be understood, however, that the 
aiming principles set forth herein may be used in connection with known TMS delivery 

30 modes, and likewise, the novel stimulator and delivery system may be used in TMS 
delivery modes without use of the aiming principles set forth herein. 

A flat, B-shaped coil, such as is most commonly used for TMS, creates an E-field 
that is oriented tangential to the scalp surface upon which it is placed. By Aiming 
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Principles 2 and 3, this E-field configuration is optimal for stimulating the banks of sulci 
in which the sulcus is oriented perpendicular to the scalp surface (e.g., in Figure 1). The 
same principles, however, predict that this E-field configuration will be unable to 
stimulate the cortical surfaces parallel to the scalp, i.e., the gyral crowns. The gyral 
5 crowns compose approximately 30% of the total cortical surface. Thus, there is a need 
for a TMS coil capable of stimulating gyral crowns. By Principles 2 and 3, this would 
require an E-field oriented perpendicular to the gyral crown, i.e., essentially 
perpendicular to the scalp surface. 

FIG. 7 is an example embodiment of a TMS stimulator capable of stimulating 

10 gyral crowns, according to present invention. As shown in FIG. 7, TMS stimulator 100 
includes a body portion 1 10 and a coil portion 120. In an example embodiment, the body 
portion 1 1 0 may be made of air, ferrite or other materials with \x > 1 . Stimulator 1 00 may 
be surrounded by a conducting fluid. The conducting fluid can fill the central portion of 
the coil portion 1 20. The conducting fluid can fill the space between the coil and targeted 

15 body tissue. The purpose of the conducting fluid is to provide efficient induction of the E 
field produced by stimulator 100 through the tissue surface, directed according to the coil 
design. In certain embodiments it may be desired to construct the body portion 1 10 of a 
non-conducting material and include an insert through at least a portion of the hollow 
portion of the cylinder which may be made of a ferromagnetic material such as iron. The 

20 coil portion 120 may be made of wire, such as copper or the like. Via conductors (not 
shown in FIG. 7), the coil portion 120 is connected to a power supply. The power supply 
is used to provide a current through the coil, which generates an electric field emanating 
from the stimulator. 

In an example embodiment, the power supply may include a high energy 

25 capacitor bank, which when discharged provides a high current through the coil portion 
120. In an example embodiment, the electrical current may be relatively strong (between 
about 1000 and 2000 amps) and last for a relatively short period of time (between 50 and 
250 microseconds). To deliver such a high intensity, short pulse current, the power 
supply may include a capacitor bank having a very high storage capacity on the order of 

30 50 microfarads charged to a potential of 2000-6000 volts. Preferably, the capacitor bank 
is made of high duty material to withstand the current generated. For example, the 
capacitor have a capacitance of 50 MFD and a voltage rating of 7.5 kV may be made 
physically larger to extend the lifetime from an estimated 10 pulses to 10 pulses. 
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Furthermore, in an exemplary embodiment, the power supply, conductors, and the TMS 
stimulator may be water-cooled to reduce operating temperature. 

Although shown in FIG. 7 as a cylinder, it is to be understood that body portion 
1 1 0 may be various shapes, such as cone shaped, cylindrical, ellipsoidal, rectangular, and 
5 the like. Further, as shown in FIG. 7, the coil portion 1 20 is preferably toroidal shaped, 
however, it is to be understood that the coil may be constructed in various shapes. 

In an example embodiment, the body portion 110 may comprise a cylinder having 
a diameter between about 10 and 20 centimeters, and a hollow core having a diameter 
between about 0.5 and 2.0 centimeters. 
10 In an example embodiment, the wire portion 120 extends substantially around the 

body portion 1 1 0 in a toroidal shape. In certain embodiments, the wire used may have a 
. diameter between about 0.1 and 1.0 millimeters. 

By use of such a design, when energized by an electric current, the TMS 
stimulator 100 generates a maximum induced electric field extending along the long axis 
15 1 30 of the body portion 1 1 0. 

In an example embodiment, TMS delivery may be effected via use of a robot. 
More specifically, a robot having six or more degrees of freedom may be used to 
appropriately position and orient the TMS coil in a precise location and orientation for 
most effective delivery of induced electric field. Without a sixth degree of freedom, the 
20 orientation aspect of the CCCAP described above cannot be implemented. In an example 
embodiment, a surgical robot having five joints may be modified to add a device to hold 
the TMS stimulator and provide for a sixth degree of freedom. Alternately, a robot 
having six degrees of freedom may be used. 

The robot may be a medical robot to provide benefits of a robot safe and effective 
25 for use with human patients. Such a robot may be controlled by a computing system, 
such as a PC. Robot tool motion and orientation may be controlled from appropriate 
software programs. Such programs may include control for robot arm position 
monitoring, position modeling, plan movement modeling, and movement model 
execution. Additionally, sensors may be provided so that XYZ position and angles for 
30 each of the joints of the robot arm may be obtained in real time. Motion commands may 
be executed in a coordinate system native to the robot. In an example embodiment, the 
robot arm may be moved from point to point in a smooth manner using a 3-D mouse, 
SPACE MOUSE™, or keyboard input. By use of such a robot, maintenance of a pre- 



.19- 



WO 02/089902 PCT/US02/14157 
specified treatment position may be more precise and accurate than a passive arm or a 
hand-held device. Furthermore, a robot allows movement through a treatment zone 
(rather than treatment of a point) and further allows precisely timed treatments of 
multiple sites in a single session. 
5 Via use of such a robot, TMS may be effective for use in connection with pre- 

operative mapping, inter-operative monitoring, and treatment of conditions, such as long- 
term depression, and for clinical activity and functional mapping using TMS and PET. 
Additionally, by use of a robot, a means is provided to export a generic treatment plan. In 
other words, a coordinate system may be developed and standardized across different 
10 patients, labs, and the like. Furthermore, such use of a robot is particularly useful with 
regard to a treatment program in which a patient is subjected to TMS over a repeated 
number of days. 

In an example embodiment, a NEUROMATE™ stereotactic assistant system 
(ISSI, Sacramento, CA) was modified for use. FIG. 8 shows an overview diagram of the 

15 NEUROMATE™ system 300 as used in connection with PET imaging. The exemplary 
surgical NEUROMATE™ robotic arm was modified to convert it from a system for 
frame-based, intra-operative probe positioning into a system for frameless, extra- 
operative, positioning and orientation of an induced E-field. 

Robotic devices must be used with caution in the presence of persons. In an 

20 example embodiment, a robot such as the NEUROMATE™ specifically designed for use 
by and with humans is desired. Benefits include that it moves slowly, such that it can 
easily be stopped prior to any collision, loss of electrical power renders it immobile in its 
current position, i.e., it does not return to a "parking" position or make any sudden 
motions. Further, such a robot may require a key and a recessed button to actuate, being 

25 used only under direct supervision. 

A TMS-coil holder was constructed for the TMS coil to adapt to the 
NEUROMATE™. In the example embodiment, the TMS tool-holder included a tool- 
rotation axis, as lacking this, changing coil orientation would require establishing a new 
arm posture, with only a very limited number of orientation angles (at a given location) 

30 being possible. In addition, the tool holder serves as a stand-off, creating a distance 
between the final limb of the robotic arm and the TMS coil. This distance decreases 
effects of the metal arm on the coil's B and E-fields. In addition, this helps keep the 
metal arm out of the PET field-of-view, thereby lowering the possibility of PET-image 
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attenuation artifacts. Finally, the tool holder permits orientation sensing and motorized, 
computer-controlled coil-orienting. It is to be understood, however, that a tool holder 
need not include a TMS-rotation axis, and that in other embodiments (such as a 6-joint 
robot), the TMS coil may be adapted directly to a robotic arm. 
5 In the example embodiment, the TMS-coil tool holder was constructed of 

fiberglass. However, in other embodiments, the tool-holder may be made from a number 
of other materials, such as acrylic, fiberglass, delrin, or the like. As shown in FIG. 9, the 
tool holder 200 includes a proximal end 210, which may be machined to interface with 
the robotic arm tool-mounting plate 3 1 0. The distal end 220 of the holder may have a rod 
10 230 mounted to the TMS coil 1 00 and allow rotation about the coil's z-axis. For the 
example embodiment, the robotic arm tool mounting plate 310 comprises the fourth and 
fifth axes of the robot, and the tool holder 200 comprises the sixth degree of freedom. In 
an example embodiment, the body of the holder 200 is approximately 8" and houses a 
plastic belt turning a more proximally placed rod to mechanically link the orientation of 
1 5 the rod 230 and the TMS to the orientation sensor 250. 

In an example embodiment, a precision rational sensor 250 may be mounted on 
the end of the proximal (belt-driven) rod. In an example embodiment, the sensor may be 
a single turn (340°, lOkQ) precision resistive potentiometer with +/- 1% linear rating. 
The output of the sensor specifies the orientation of the TMS coil. The sensor's output 
20 may be conveyed to TMS delivery software via tool-sensing electrical contacts on the 
robotic arm tool-mounting plate 310. In neurosurgical applications, these contacts are 
used to monitor the depth of tool (e.g., drill-bit or wire electrode) passed through a 
stationary tool holder. In the present application, they may be used to monitor orientation 
rather than depth. Additionally, a stepping motor may be added, allowing active rotations 
25 of the TMS coil. 

As discussed above, the NEUROMATE™ was designed as a system for framed 
stereotaxy. A patient is registered to the NEUROMATE™ by having his (or her) head 
placed within a metal ring that is rigidly mounted to the NEUROMATE™ by means of a 
pedestal attached to a long leg, which forms the base/stand of the NEUROMATE™. For 
30 TMS delivery, the robotic arm was converted into a frameless stereotactic system, 

capable of being rapidly registered to a subject in any position within the operating range 
of the robot arm, including supine in the PET scanner. It is to be understood, however, 
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that the methods and apparatus of the present invention may be used in conjunction with 
framed stereotaxy. 

In an example embodiment, this registration process may use a manual digitizer to 
collect points from the subject's scalp, which may be registered to the scalp surface 
5 modeled from that subject's image, such as obtained from MRI. In other embodiments, 
- the robot may be used in an interactive mode to collect these scalp points. 

In an example embodiment, a 5-axis manual digitizer (Microscribe 3D) may be 
permanently mounted to the robot chassis for use as required. The positional relationship 
between the digitizer and the TMS probe may be calibrated and thereafter should remain 
10 constant, although it is recommended that calibration be verified for each use of the 
system. 

Typical surgical robotic devices come equipped with a stand designed for framed 
stereotaxy. This stand has a long foot which prevents it from being placed sufficiently 
close to the PET for use during TMS/PET. Thus, a shallow, wide, weighted, four- 

1 5 wheeled, cart may be used having weight and width to prevent tipping. Preferably, the 
wheels are retractable, to allow the system to be completely stable when in operation, yet 
re-positionable to optimize placement of the robotic arm relative to the PET. 

As presented above, passive digitizing arms can be used for image-guided TMS- 
delivery. Both Paus (et al., 1997) and Krings (et al., 1997) used commercially available 

20 passive arms designed for frameless stereotaxy to assist with TMS aiming. However, 
these two systems suffer from significant and virtually identical limitations. Neither can 
be fixed rigidly in position; that is, both are manual digitizers but neither is a holding 
device. (Paus mounted the TMS coil to a separate, home-made holding arm; Krings held 
the TMS coil by hand.) Neither system can digitize or visualize tool orientation, as both 

25 are five-joint arms with no capability for rotation about a tool axis. (Paus fixed the coil in 
the B-nose orientation; Krings hand-held the coil in the orientation recommended by 
Brasil-Neto [et al., 1992] and Mills [et al., 1992]). Not unexpectedly, neither tool has 
software for visualizing tool position/orientation vectors. Neither system has any 
capabilities for modeling the E-field, visualizing E-field lines relative to brain anatomy. 

30 Finally, neither system is robotic. 

An active (robotic) arm has clear long-term advantages. Ideally, treatment 
positions are pre-specified, with the arm moving to and maintaining the treatment 
position. A robotic system can carry out this function more precisely and more 
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accurately than a passive arm. Further, only a robotic system readily combines the tool 
aiming and tool holding functions. Finally, only a robotic arm allows movement through 
a treatment zone (rather than treating at a point) or allow precisely timed treatments of 
multiple sites in a single session. 
5 The PET image, the MR image, the manual digitizer, and the robotic arm robot 

have different coordinate-reference frames. In an example embodiment, a whole-brain, 
high-resolution (1 mm 3 ), Tl -weighted MR! may be used to create a reference space 
within which all objects become registered. This image format provides ample 
anatomical detail for precise registration. Additionally, all images may be standardized 

1 0 by removing differences in orientation and position between patients. In an example 

embodiment, each patient's reference MRI may be rotated and translated (but not scaled) 
to the alignment of the Talairach and Tournoux (1988) atlas. This six-parameter spatial 
normalization (of the reference MRI) may be performed using the previously validated 
SN software (Lancaster, et. al. 1995). This software may be used to freely rotate and 

15 translate a 3-D brain image to position it in a standard pose. The standard pose is with 
the brain's anterior commissure (AC) located at x, y, z, — 0,0,0. The line between the 
anterior and posterior commissures (AC and PC) adjusted to lie along the y-axis and the 
inter-hemispheric plane fitted to the y-z plane at x = 0. The reference MRI may serve as 
the reference frame for coordinates from: functional images (PET), the head-surface 

20 digitizer; the robotic aim, the TMS coil (held by the arm), the standardized head models; 
and the TMS-induced B and E fields. 

In an example embodiment, the robot may be moved in a user-guided visual- 
feedback mode for manual positioning and holding of the TMS probe. Preferably, the 
robot may be moved from point to point in a smooth manner using joy stick-like or 

25 keyboard input. For example 3-D mouse/joystick, such as a SPACE MOUSE™ may be 
used in conjunction with a supervisory PC controller to send motion commands to the 
robot's on-board controller. Motion commands may be executed relative to: (1) the tool 
coordinate system referenced to the face of the TMS; (2) a reference MRI (aligned but 
not scaled to the Talairach & Tournoux [1988] brain); and (3) the global coordinate 

30 system native to the robot. In order to implement this control feature, the supervisory 
controller may use a robot inverse transformation to calculate joint coordinates needed to 
accomplish an incremental change of position relative to tool coordinates. The 
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calculations may take into account arm configuration and joint rotation limits, and will 
format the commands to the world space required by the supervisory controller. 

In an example embodiment, the communication protocol, position-sensing 
interface specifications, movement sensing interface specifications and complete 
5 kinematic equations for the robot arm may be developed using, for example, Microsoft 
Visual C++ on a Pentium-level PC. 

The protocol used by the supervisory control software to communicate with the 
robot controller may use RS-232 connections and protocol to send instructions and 
retrieve supervisor-level data from the controller. Using such protocols, a Visual C++ 
10 application may support all motion type commands, including an exclusive stop 
command and polling commands for status. 

A second serial port of the PC may be used to get robot arm position data and 
status. The communication protocol may be similar to that for the supervisory control 
software discussed above. A multi-port serial I/O board may be added to the PC to 
1 5 support simultaneous access to both ports from the same computer. Supervisory control 
software routines may be used to read this port and obtain x-y-z position and angles for 
each of the six joints of the robot arm. Furthermore, forward and inverse kinematic 
equations may be used to specify and modify arm pose. 

In an example embodiment, the patient's head and the robot may be registered to 
20 the MR reference frame for each patient and each TMS session, as accurate positioning of 
the robot forms the basis for registration of the TMS coil relative to the head. Physical 
coordinates digitized about the surface of the head may be used to register the patient's 
head and to calibrate the digitizer relative to the MRI reference frame. As disclosed 
above, the manual digitizer may be permanently mounted on the robot, enabling direct 
25 digitizer-to-robot calibration. Rigid-body coordinate-transformation matrices may be 
calculated between the robot, the digitizer, and the MR reference frame. 

For patient registration, a manual 3-D digitizer (MicroScribe 3DLX) may be used 
to acquire a "cloud" of distributed points on the patients head in a hat-file format 
(Pelizzari et al., 1989). An MRI surface (head-file format) with 6000-7000 points 
30 uniformly distributed about the head may be extracted from the subject's standardized 
MR image using, for example, convex hull (CHSN) registration software (Lancaster et al. 
1999). A coordinate transformation for calibration (or registration) may be determined 
by fitting the digitized hat "cloud" (patient) to the detailed head (MRI) surface. The 
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fitting method used may be an iterative least square technique that has been shown to 
provide excellent results for co-registering PET, CT, and MRI (Lancaster et al., 1997). 
Patients may be positioned within the PET scanner with head constrained. The CHSN 
software may be used for fitting the hat to the head surfaces and calculating a head-to- 
5 MR transform matrix. The digitizer-to-MR matrix is identical to the head-to-MR matrix, 
as both are calculated using the manual digitizer and the head surface. 

For robot registration, the robot-to-digitizer transform matrix may be calculated 
using fiducials for which both digitizer and robot coordinates can be readily determined. 
A pseudoinverse method may be used to calculate the transformation matrix (Castleman, 

1 0 1 996). The robot-to-MR transform matrix may be calculated by concatenating the 
transform matrices for robot-to-digitizer and digitizer-to-MR. 

When the TMS coil is permanently mounted on the robot arm, the TMS-to-robot 
calibration transform matrix can be calculated. This may be done with the coil in a 
designated zero-degree orientation relative to the mount. Three parameters (x-y-z 

15 position, coil, tool, or z axis orientation, and rotation about coil, tool, or z-axis) are used 
to fully aim the TMS coil. These may be determined using the manual digitizer to record 
fiducials on the coil's x- and y-axes. Using these data, equations to calibrate the robot 
can be determined and robotic coil aiming enabled. This calibration should remain fixed 
but may be verified periodically. Coil aiming may be extended to the MR reference 

20 frame using the robot-to-MR transform matrix. Only position and z-axis orientation are 
under direct control of the robot. The uncontrolled rotation about the z-axis may be 
calculated for each robot arm position. Targeted orientations of the coil about its z-axis 
may be achieved using the manual digitizer and/or the orientation sensor. 

In order to overlay PET (functional) and MRI (anatomical) images as part of 

25 treatment planning, the PET images may be registered to the MR reference frame for 
each patient. PET images may be co-registered to standardized MR images using, for 
example, the convex hull (CHSN) registration software (Lancaster et al., 1999). Convex 
hull surface is extracted from both MRI and PET images. These models remove concave 
regions that are well resolved by MRI but poorly resolved by PET. However, the convex 

30 brain surface is identical in both imaging modalities and proves and excellent surface for 
surface matching and, thereby, volume registration. The CHSN software (Lancaster et 
al., 1999) uses the convex hull of the brain, a simplified surface which is highly similar 
across imaging modalities, for registration and spatial normalization. An iterative least 
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square method is used to fit the convex hull from PET (hat format) to the subject's 
reference MRI (head format), allowing only rotations and translations of the hat data 
(Pelizzari et. al., 1989). This fitting method has very small error: mean RMS errors, 
measured as distance between the two convex hulls, were reported to be less than 1 mm 
5 (Lancaster et al., 1 999). A PET-to-MR coordinate transformation matrix may thus be 
calculated by the CHSN software. 

There are four principles that can be used to estimate important aiming parameters 
of the TMS stimulator without resorting to comprehensive E-field calculations. First, the 
E field magnitude of figure-8 coil tends to remain strongest at points along the central 

10 axis perpendicular to the plane of the coil (z-axis). Second, the E-field magnitude rapidly 
decreases with increasing distance from the surface of the coil. Third, inside a conductor 
(e.g. brain) the surface charges reduce the E-field magnitude relative to what its 
magnitude would be at the same point in air. Four, under certain conditions of symmetry 
a weak conductor (e.g. the brain) will not change the direction of the E-field on the coil's 

1 5 z-axis. For example, if the coil is placed tangential to the head, at a point of fairly 

symmetric curvature, this symmetry about the z-axis will tend to cancel distortions that 
would change the direction of the E-field. These "on-axis" principles, in conjunction 
with the columnar aiming principles may be used to accurately aim the TMS. 

To create a safe and effective treatment plan, the following exemplary steps may 

20 be performed. First, a cortical stimulus-delivery site must be selected. Second, at the 
stimulus-delivery site, the columnar orientation must be determined. Third, TMS coil 
aiming parameters must be determined. Fourth, a model of the TMS E-field (corrected 
for head-induced distortions) must be determined. Fifth, models of the TMS coil and 
patient head surface must be provided for collision avoidance. However, it is to be 

25 understood that in other embodiments, more or fewer steps may be performed. 

In an example embodiment, these steps may be performed via TMS planning 
software. Site selection may be done interactively while viewing a high-resolution MR 
image of the brain in three orthogonal views. The columnar orientation at the treatment 
site may be determined from the high-resolution MRI by means of cortical-surface 

30 extraction, as cortical columns are oriented perpendicular to the cortical surface. Once 
the treatment site and its columnar orientation are established, the on-axis plan proceeds 
by determining the TMS coil position for aiming. This may make use of a 3-D front- 
surface model of the coil and the patient's 3-D head surface model to avoid head contact. 
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The coil position closest to the head in which the on-axis E field parallels the cortical 
columns at the treatment site may be established. The TMS coil aiming parameters 
(position, coil z-axis orientation, coil y-axis orientation) may be determined and used to 
direct the robot to aim the coil during TMS delivery. 
5 Using the TMS coil aiming parameters and the treatment site configuration 

(location, column orientation), the induced E field may be modeled using on-axis data to 
predict its magnitude at the delivery site. A fully 3-D approach, including effects of the 
patient head may be modeled for each subject using EM AS (Ansoft, Pittsburgh, PA) and 
anatomical data derived from CT scans. 

10 In still further embodiments, methods for accelerating calculations of electric 

fields may be provided. EMAS is a general purpose Maxwell's equations solver. 
Because of its generality, it is slow and uses a lot of computing resources. Calculation of 
the electric field for a single coil position may require several hours for interactive, 
image-guided imaging abbreviation of this computation time is desirable. Computations 

1 5 may be accelerated chiefly by introduction of simplifying assumptions. For example, 
because the skull is a poor conductor, the magnetic fields produced by induced eddy 
currents can be safely ignored (Davey et al, 1991). As a consequence, the magnetic field 
produced by the TMS coil may be modeled as being unaffected by the head. Further, the 
skull is 100 times less conductive than soft-tissue. Therefore, it may be possible to treat 

20 the skull as a perfect insulator; treating all regions outside the skull as air. In addition, 
when using a standardized head model, it may be possible to pre-compute (and store as 
look-up tables) the E-field for often-stimulated brain zones, such as Ml -hand. The E- 
fields for an individual could be then be calculated as a perturbation of the pre-computed 
values. Finally, it is possible that field estimates based just on the local skull curvature 

25 may be quite accurate and provide considerable computational savings. 

To support safe movement of the coil about the head, 3-D models of the head and 
TMS coil surfaces may be defined, and a polygonal model of the surface of the TMS coil 
created. Full-surface polygonal models of each subject's head may be created from the 
high-resolution 3-D MR image. By thresholding the MR image to the head surface, a 

30 mask defining the full 3-D volume of the subjects head can be obtained. Using 3-D 
morphological dilation, a second volume that is guaranteed to be a minimum distance 
from the head surface can be created. The dilated head volume may be used to model a 
volume with a realistic safety margin for use in collision avoidance control when moving 
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from point to point about the head. The undilated model may serve to determine closest 
approach without contact in final positioning. The surfaces of the head models may be 
extracted and stored using the marching cubes algorithm and its standard storage format 
(Schroeder et. al., 1996). Additionally or alternatively, collision avoidance may be 
5 provided by the operator (run/stop button with absolute shut-off switch). 

In an example embodiment, TMS delivery software may be used to position the 
TMS coil exactly as prescribed by the planning software. This software may provide 
accurate robot-controlled positioning of the TMS coil along and about its z-axis. 
Additionally, movement and pose strategy may ensure safe, effective operation. Finally, 

10 feedback to the delivery software may monitor this operation. 

In exemplary software, a graphical user interface may be provided, in which the 
user is prompted for input and presented with options for delivery of TMS. The TMS 
delivery may proceed and display status information. A graphical user interface may 
monitor movement of the robot/TMS coil about the patient's'head. Polygonal models 

15 developed for collision avoidance may be used for rendering 3-D views of the head and 
TMS coil body. A 3-D polygonal model of the robot with TMS coil mounted may be 
used to create a real-time 3-D display of the robot moving about the head surface. View 
orientation may be user-selectable to provide unobstructed viewing of the patient's head. 
This display feature may be provided as animation during treatment planning to simulate 

20 delivery. 

For a given TMS coil location and orientation, several arm poses will often be 
possible. Environmental factors, including TMS coil electrical and cooling lines, PET 
scanner housing, patient table and head holder, etc., may be accommodated when 
selecting best pose. Joint configuration may be under control of the supervisor software 

25 developed to move the arm, which are computed using the forward and inverse 

kinematics equations of the robot. The user may be presented with several options and 
asked to select one for the treatment delivery. 

In an example embodiment, the robot can position the TMS coil but cannot rotate 
it about coil's z-axis, therefore the coil angle may be adjusted manually. Software may 

30 use the manual digitizer to measure the TMS coil z-axis rotation angle. In embodiments 
with a rotational sensor, the sensor may be used for coil-angle sensing, in preference to 
the manual digitizer. The current and target angle, and the difference may be actively 
displayed. Manual orientation may be done with TMS housing away from the head at the 
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position where the robot has already aligned the coil along its z-axis. The z-axis rotation 
angle sensing may be automated using the precision rotary potentiometer. The sensor 
may be connected to the robot controller via its existing wiring harness. Movement 
toward the head along this treatment axis may be done after adjusting coil angle. 
5 The cortical column orientation throughout the treatment volume may be 

determined by a method using both inner and outer cortical surfaces. The first step may 
be a semi-automated extraction of the cortical surface, using, for example, the MEDx 
(Version 2.1, Sensor Systems, Inc.) image processing software. This step requires 
approximately 20 minutes per surface. The extracted surface may be refined using 

10 morphologic operators (dilation, erosion, and 3-D connectivity routines), also within 
MEDx. Both inner and outer cortical surfaces for the treatment volume, including all 
visible sulci, may be extracted and used to create a binary mask of the cortical columns 
within the treatment volume. The binary mask may be used to create a set of normal 
vectors that uniformly fill each voxel within the mask. A 3-D Euclidean Distance Map 

1 5 (EDM) (Russ 1 996, Gonzales and Wood 1 992) may be used to fill the binary mask with 
distances from the outer to the inner surface. The EDM may then be smoothed with a 3- 
D Gaussian filter to produce a smooth gradient of values from the outer to inner cortical 
surface bounds. A normal vector may then be calculated for each voxel within the binary 
mask from the gradient of EDM values (Schroeder, et. al. 1996). These normal vectors 

20 may be converted to unit normal vectors. The sense of the unit normal vectors may be to 
point into the brain (orthodromic columirdirection). This method, based on the EDM 
provides a smooth and continuous orientation change for the normal vectors. These 
normal vectors may be used for the orthodromic orientation of the cortical columns for 
each voxel in the binary mask. 

25 A 3-D scalar map of columnar component of E-field may be accomplished using 

3-D models of the E field accounting for the head a set of coil aiming parameters 
(position, coil z-axis orientation, and coil rotation), and the 3-D model of cortical column 
orientations. Using the TMS coil's treatment parameters to orient the E-fields relative to 
the head, the scalar product of E-field vectors and cortical-column unit vectors can be 

30 calculated, yielding E-field magnitude along cortical columns (Volts/cm) for each voxel 
within the treatment volume. The result will be a 3-D scalar map of the columnar 
component of the E-field throughout the treatment volume. This map may be used to 
overlay onto the MR image to see the planned treatment. Modeled treatment effects may 
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be visualized as a pseudo-color overlay on the anatomical image. An option may be 

provided to display the absolute value of the E- field magnitude, and a method to 

enable/disable the overlays. The user may adjust TMS coil global E-field magnitude and 

a threshold value for display to see the effect of these on target and non-target cortical 

5 regions. Since the E-field scalar map depends on both cortical column and coil 

orientation, a "jog" feature to interactively reposition the TMS coil may be provided to 

help the user optimize the magnitude in a targeted region. This jog feature allows small 

increments of movement along and rotations about the orthogonal tool axes after they 

have been verified safe (non-interfering with the avoidance volume) by the supervisory 

10 controller. 

In certain embodiments, the metal in the PET gantry might alter the TMS-induced 
B and E fields. If so, this effect must be included in the E-field modeling. Conversely, 
the TMS B-field might affect the performance of the PET photo-multiplier tubes, which 
may be resolved with \i metal shielding during TMS/PET. Further, the TMS coil and coil 
15 holder may excessively, regionally attenuate the emitted annihilation photons. It is to be 
understood that standard attenuation corrections may suffice to correct for the effect of 
the TMS coil. 

FIG. 10 is a block diagram of an example TMS treatment planning system 
according to the present invention. As shown in FIG. 10, in an example treatment 

20 planning system 400, a desired site is defined at step 410. For example, the desired site 
may be defined via various imaging techniques, such as fMRI or PET. Next at step 420, 
the TMS delivery parameters may be determined. These parameters may include, for 
example, the desired pulse rate, duration, TMS stimulator design, and a definition of 
critical/non-critical regions. Cortical regions may include sites to avoid during treatments 

25 or regions of known connectivity to targeted regions that might alter the net stimulation to 
the targeted region. Non-critical regions might be regions known to have little 
connectivity to the targeted region. 

At step 430, the defined site may be registered to the subject's anatomy. In an 
example embodiment, the registration step may be performed via use of a 3-D MRI. At 

30 step 440, planning tasks may be performed. These planning tasks may include modeling 
various elements in order to optimize TMS delivery. These models may include, for 
example, a 3-D cortical model of the subject's brain, to determine the accurate location of 
cortical columns and cortex regions. In addition, a 3-D E-field may be modeled for the 
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chosen TMS stimulator design. Further, in an example embodiment, a 3-D induction 

model may be performed also. The 3-D induction model is derived from the scalar 

product of the E-field and the cortical column direction, and is the net volts/cm estimated 

for cortical columns. Based on all this information, at step 450, the TMS treatment plan 

5 may be computed. In an example embodiment, the plan may include information 

regarding optimal coil position, coil orientation, power setting, rate/duration, and a 3-D 

head model. The TMS delivery plan developed at step 450 may then be provided to a 

TMS treatment delivery system. It is to be understood that the TMS treatment planning 

system described herein may be performed using a standard personal computer 

10 appropriately programmed, or it may be performed via specialized computers, such as a 
UNIX workstation or a mainframe computer. 

FIG. 1 1 is an example flow diagram of a TMS treatment delivery system 
according to the present invention. As shown in FIG. 1 1 , an example TMS treatment 
delivery system 500 may begin at step 510, in which a patent is positioned with his head 

1 5 immobilized so that a duplicate orientation (the same orientation as used for imaging) 
may be obtained. At step 520, spatial calibrations may be performed utilizing the data 
provided by the treatment planning system. Spatial calibrations may include, for 
example, a robot to digitizer calibration, a digitizer to MRI calibration, and a robot to 
MRI calibration. In an example embodiment, the robot digitizer spatial calibration may 

20 be performed at step 530 in which a physical coil model is acquired via a digitizer and is 
matched with the robot coil model. Key landmarks on the coil body may be used to 
establish coil x-y-z orientation, the coils' z-axis direction, and the 0,0,0 point on the coil 
surface. At step 540, the digitizer to MRI calibration may be performed in which a 
physical head model is acquired via the digitizer and is matched with the 3-D MRI head 

25 model provided by the TMS treatment planning system, as discussed above. The robot- 
to-digitizer and robot-to-MRI transforms may be concatenated to create a robot-to-MRI 
transform during calibration. Next, at step 550, safety features may be implemented in 
which a plan for avoidance of the head during movement of the robot is effected. 
Furthermore, the rate/duration provided by the TMS treatment planning system may be 

30 implemented as a control to prevent excess stimulation. Finally, at step 560, TMS 
delivery is effected after the coil is positioned and oriented and the rate, duration, and 
power of the delivery is set. The TMS treatment delivery system may be performed using 
the computing systems discussed above. 
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In conjunction with the CCCAP disclosed herein, a higher degree of control over 
the electrical field induced by a TMS coil will improve the value of the 
experiment/treatment. Such control may be aided by new and inventive coil designs. 

Greater control may be had over such features as coil inductance and power 
dissipation with methods that have been in use for designing gradient coils for magnetic 
resonance imaging ("MRI") since the late 1980's for cylindrical coils (Turner) and later 
for planar coils (Martens et ah). These methods involve an inverse technique which 
allows the magnetic field to be defined at certain points in a volume. This method may 
be adapted to the design of magnetic stimulation coils which enables the electric field to 
be defined at certain points. While it is not reasonable to focus the magnetic field at an 
arbitrary position inside the brain/muscle volume, the resultant coils may provide more 
focused fields in two dimensions near the brain/muscle surface and with other desirable 
characteristics such as a lower inductance or a lower heat dissipation. 

Suggestions of 10 kHz as the optimum frequency for the stimulation waveform 
have been made (Davey and Epstein). The stimulation circuit may be considered 
primarily an LC circuit with a small resistive loss. The capacitive section of the circuit 
stores energy and when a switch is closed, the energy is transferred to the inductive 
section, the magnetic stimulation coil, and back to the capacitive section with a frequency 
<D=l/(LC) l/2 . The voltage induced in a conducting body is proportional to this frequency. 
For the purposes of coil design, the coil is of interest, and thus the capacitive section may 
be considered constant. Hence, to increase the resonant frequency and the induced 
voltage, the coil inductance may be lowered. 

The target field method may be adapted to provide a coil with a minimum 
inductance for a given set of electric field constraints. The constraints may be set to 
provide as smooth a design as possible. The constraints should also be set such the coil is 
of a reasonable size in order that the natural field falloff moving away from the coil is not 
too sharp. 

The theory behind minimum inductance TMS coils may be adapted from known 
design techniques for target field gradient coil designs for MRI, such as that of Martens et 
al, and Turner. For current flowing on a planar surface, the x-y plane, the Fourier 
transform pair may be defined as: 




(i.i) 
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and 



Jx (*> ') = jj e""e'»j, (a, p)dadp (2) 



and similarly for the y-component of the current density. Using: 
div J = 0 (3) 
5 leads to, at the plane y = a: 

cg„(a,p)+ff M (a,t)=0 (4) 

From the vector potential: 

4n JJv |r-r | 
Orthogonal components may be determined: 

and 

4(r) = 0. (8) 
The electric field in a homogeneous media may be written as: 

15 E = -^ (9) 
ot 

and so for a sinusoidal current, I = 1 0 sin o)t, 
E x (r) = A x (r) a> cos cot etc. (10) 

20 If the x-component of the E-field is considered: 

E x\=®- A X (11) 
and so 



(6) 



(7) 
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K\ =°> A ] ^-«*«- P5r 'l/,(«.^|rf«* 

^ -co — oo'Y^ 

(12) 

The electric field must be symmetric in both the x and y directions, so that the imaginary 
components in the above equation can be reduced to cos(ax) and cos(Py). The energy 
from the current density plane may be calculated from: 

W = i Ja.JJ 3 * (13) 



2v 



leading to: 



16* 2 i} Ja 2 +/3 2 { 

In order to minimize W subject to a set of electric field constraints, the error function may 
10 be formed: 

U[j;(a,p)] = W-X^j[ E ( r j)- Er ( r i)] O 5 ) 

H 

where \ } are Lagrange multipliers and the electric field at points rj is constrained to have 
the values E l (r,). Setting: 

dv 



5j x a (a,P) 
15 gives: 



= 0 (16) 



2 

(17) 



a + p jfi 

The Lagrange multipliers are found by solving the matrix equation: 

E'Cr^JCyXj fori = l,N (18) 
j=i 

found by substituting the current density equation into the electric field equation. 

20 Once the current density has been calculated, the Fourier transform may be used to obtain 
the current density J x (x,z) and J 2 (x,z). These may then be integrated to provide a stream 
function and may be discretized into current loops. The inductance may be estimated 
using W = 0.5*L.I 2 . The accuracy of the design is then dependent on the number of loops 
chosen to approximate the continuous current density. The design can then be checked 

25 using: 
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(19) 



and summing individual wire elements to obtain the total field distribution. The coil may 
also be evaluated in the presence of conducting boundaries to simulate a human head 
using numerical methods such as finite element analysis. The design process does not 
5 take into account these boundaries, although an approximation may be incorporated into 
the design using a conducting sphere and a well-known method of calculating electric 
field in a spherical volume, as set forth in Eaton (1992). 

To obtain a varying electric field, the symmetries may be adjusted, i.e. all cos(ax) need to 
1 0 be changed to sin(ax). 

To achieve a minimum power design, the above calculation may be used, but replacing 
the energy term with the following: 



15 potency of magnetic stimulation. For situations involving rapid stimulation, heating of 
the coil becomes an issue due to the increased duty cycle. While cooling mechanisms 
have been employed it is also advantageous to include in the coil design a power 
dissipation term, in order for heating to be minimized. Minimum power designs are 
larger and generally smoother than minimum inductance designs. 

20 When attempting TMS, unwanted stimulation of the nerves in the scalp may occur 

for persons with high thresholds (usually meaning a slightly thicker skull), possibly 
causing pain and facial twitching, confounding experimental results. With the addition of 
negative turns to the windings of the stimulator coil, and a small spacing between the coil 
and the scalp, the fields across this area may be reduced significantly. However, this 

25 comes at the expense of some focusing ability, and coil inductance and heating. 

FIG. 12 shows an example embodiment of a coil design developed by an inverse 
method. As shown in FIG. 12, coil 600 is a two-winged coil design having a first wing 
620 and a second wing 630. The wings extend peripherally from a center portion 610. 
As can be seen in FIG. 12, coil 600 is a variation of a figure-eight coil design. The coil 

30 600 varies from such a standard design, in that the individual wire elements or windings 




(20) 



Recent studies have suggested that rapid stimulation (> 1 Hz) may enhance the 
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640 that make up the coil are closely concentrated towards the center portion 610 and 

extend therefrom in increasing arcs to the periphery of the coil. As shown in FIG. 12, the 

direction of the electric field traveling through the windings 640 is opposite directions for 

each of the two wings of the coil. However, the windings 640 making up each of the 

5 wings 620 and 630 all travel in the same direction. 

The spacing between windings 640 may vary in given embodiments. As shown in 

FIG. 12, there are approximately 12 wires in each wing of the coil. However, it is to be 

understood that more or fewer windings may be present in a given design. Further, the 

wings 620 and 630 are generally shown as mirror opposites. However, it is to be 

10 understood that in certain embodiments, the wings may have different eccentricities. 
The coil may be designed in accordance with standard practices for coil 
manufacture well-known to those with skill in the art. In an example embodiment, the 
windings may comprise copper or another conducting material, such as silver (or other 
conducting material) ribbon placed on its edge. Such wire may have a diameter of 

15 between about 0. 1 and 1 .0 millimeters in example embodiments. Furthermore, the coil 
may be encased in plastic or another non-conducting material in order for reasons of 
safety and other issues. In certain embodiments, the coil may be encased in thermally 
conductive epoxy to enhance heat dissipation. Further, the coil may be encased with a 
water, oil or air cooling system. Coils may also be wound in layers placed directly upon 

20 one another, and connected in series or parallel. The wire pattern may be transferred to a 
sheet former for construction using a computerized milling machine, or may be 
transferred using heat from an inked hard copy, such as a computer printout. The wire 
pattern may be constructed from one single wire/ribbon wound with small connecting 
paths between loops. Standard connections from the coil to cabling necessary to adapt 

25 the coil to a magnetic stimulator may be made. 

The method may be changed slightly to provide a minimum inductance 
configuration providing a spatial gradient in the electric field for peripheral nerve 
stimulation. An example of this is shown in FIG. 13. 

FIG. 13 is a top view of an alternate coil design according to the present 

30 invention. As shown in FIG. 13, coil 700 is a four- winged coil having four wings 720. 
Each of the wings 720 has a first set of windings 722 and a second set of windings 724. 
As seen in FIG. 13, each of the wings 720 is generally a variation of a figure-eight coil, in 
which the more centrally located portion 722 is smaller than the peripherally located 
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portion 724. This coil may be designed similarly to coil 600 discussed above, in that each 
of the portions of each wing 720 may be designed so that the windings are concentrated 
more closely towards the interior portion of the coil and extending therefrom in 
increasing arcs. This is especially the case for the peripheral portion 724. 
5 In certain embodiments, it may be desirable for the coil to be designed in a 

generally square manner. The coil shown at FIG. 13 may be constructed in accordance 
with well-known principles for coil manufacture, as discussed above. 

FIG. 14 is a diagram of another coil embodiment created by an inverse design 
method. As shown in FIG. 14, coil 800 is a 6-winged coil having two central windings 

10 810 and four outer windings 815 in which the outer windings 815 extend peripherally 
from center portion 820 in increasing arcs. This coil was designed using constraints on 
the x-component of the electric field. In an example embodiment, the coil dimensions 
may be approximately 36 cm (x-direction) by 28 cm (y-direction), substantially larger 
than the B-shaped coil. However, it should be noted that in other embodiments, the 

1 5 current density may be apodized (basically, outside a certain region the current density 
may be multiplied by an exponential function so that the outer rings are effectively 
shrunk). The electric field may then be recalculated to determine the significance of the 
effect. 

The constraints in this embodiment were set such that the x-component of the 
20 electric field at the center point was twice that of a point at 2.5 cm along the x and y axes 
and 3 cm from the coil. That is, the constraints were: 
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FIGS. 15 and 16 are graphical representations of the electric field magnitude/current 
along the x and y axes, respectively. The current required to obtain a field of 100 V/m at 
the central point was 2951 amps, assuming an excitation waveform frequency of 5 Khz. 
The coil inductance was estimated to be 20uH. The three constraints (one at the center, 
30 one each on + x and y axes) represent an attempt to localize the generated field inside a 
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given region. Other constraints (- x and y axes) are implied by symmetry. The four outer 
windings 815 serve to limit the extent of the field in the x-direction. The central 
windings 810 provide the main field generation and limit the field extent in the y- 
direction. Further, the smooth spreading of the wire paths serves to reduce mutual 
5 inductances between wires. 

FIG. 1 7 is a diagram of an alternate coil embodiment. As shown in FIG. 17, coil 
900 is a 6-winged coil having two central coils 910 and four outer coils 920. The coil of 
FIG. 17 differs from the coil of FIG. 14 in that the central coils 910 include several 
negative turns in order to reduce the magnitude of the electric field. The negative turns 

1 0 are shown more specifically in FIG. 1 8, which is a close up of the central coils 910. As 
shown in FIG. 18, the central coils 910 include a plurality of first turns 930 and a 
plurality of second turns 940. As shown in FIG. 18, the second turns 940 are configured 
in the opposite direction and carry current having the opposite polarity of the current 
carried by the first turns 930. Although shown in FIG. 1 8 as having three negative turns, 

15 it is to be understood that in other embodiments more or fewer negative turns may be 
present in a given coil design. The constraints were the same as the coil for FIG. 14, 
except for an extra constraint at x=y=z=0.5 cm, to reduce the magnitude of the electric 
field at a distance of 0.5 cm from the coil. This extra constraint resulted in a coil that 
reduced the ratio of the maximum value of the calculated electric field in a plane 5 mm 

20 from the coil to that in the desired stimulation area (3 cm away from the coil at the center) 
by 1 8%. The cun-ent needed to create the same field of 1 00 V/m is 3425 A as opposed to 
2951 for the coil of FIG. 14. The coil focusing properties are similar to that of FIG. 14. 

It is to be understood that in certain embodiments, the negative turns may be 
located in a separate layer from turns having an opposite polarity. For example, small 

25 secondary coil, consisting of a few turns of either a figure 8 design, or a smaller version 
of the target field designs disclosed herein, may be placed next to the magnetic 
stimulation coil such that the center of the coils are coincident but such that the distance 
between them is adjustable. The current may flow in the second coil such that the electric 
field generated by the current pulse is opposing the direction of the field produced by the 

30 main coil. Further, the angle that the loops make to the surface of the coil may also be 
adjusted. 
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It is to be understood that the coil designs disclosed herein may be used 

independently from the delivery methods herein. Further, it is to be understood that the 

delivery methods disclosed may be used independently from the coils herein. 

While particular embodiments of the present invention have been shown and 

5 described, it will be obvious to those skilled in the art that, based upon the teachings 

herein, changes and modifications may be made without departing from this invention 

and its broader aspects and, therefore, the appended claims are to encompass within their 

scope all such changes and modifications as are within the true spirit and scope of this 

invention. Furthermore, it is to be understood that the invention is solely defined by the 

10 appended claims. 
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WHAT IS CLAIMED IS: 

1 . A coil for transcranial magnetic stimulation, comprising: 

a solid former having a hollow portion extending therethrough on a long axis of 
said former; and 

a wire winding wrapped around said long axis and extending at least substantially 
around said solid former. 

2. The coil of claim 1, wherein said solid former comprises a field-enhancing 
material. 

3. The coil of claim 2, wherein said field-enhancing material comprises 

ferrite. 

4. The coil of claim 1, wherein said hollow portion is filled by a field- 
enhancing insert. 

5. The coil of claim 4, wherein said field-enhancing insert comprises iron. 

6. The coil of claim 1, wherein said wire winding is wrapped in a helical 
manner through the hollow portion of said former. 

7. The method of claim 1 , wherein said solid former comprises a cylinder. 

8. The coil of claim 1, wherein said solid former comprises a hollow 
truncated cone. 

9. The coil of claim 1 , wherein said coil when energized creates an electric 
field having a maximum along said long axis. 

10. The coil at claim 1, further comprising a protective layer surrounding said 
solid former and said wire winding. 

11. A TMS stimulator, comprising: 

a core element having a first end and a second end; 
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a conductor having a plurality of wrappings extending around at least a portion of 
said first end for permitting a flow of electric current for creation of an 
electric field, wherein said stimulator when energized develops an induced 
electric field having a maximum extending substantially normal to said 
5 first end; 

said TMS stimulator having a shape such that said first end may be positioned in 
close relation to a region of interest of a patient for delivery of said 
induced electric field substantially normal to said region of interest. 

1 2. The TMS stimulator of claim 1 1 , wherein said region of interest comprises 

10 cortex. 

13. The TMS stimulator of claim 1 1 , wherein said induced electric field is 
delivered substantially perpendicular to the cortical surface at a gyrus. 

14. A system for transcranial magnetic stimulation, comprising: 

a robotic member having a distal portion and a proximal portion, said robotic 
1 5 member having at least six degrees of freedom; and 

a coil for generating an electric field, said coil coupled to said distal portion of 
said robotic member. 



1 5. The system of claim 14, further comprising a computer to control 
movement of said robotic member to position said coil at a preselected position relative to 

20 a subject. 

1 6. The system of claim 14, further comprising a power supply to supply 
energy to said coil to generate said electric field, said power supply comprising a 
capacitor bank. 

1 7. The system of claim 14, further comprising a digitizer adapted to said 
25 robotic member. 

1 8. A method of transcranial magnetic stimulation, comprising: 
defining a cerebral cortical region of interest for a subject; 
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determining a stimulator position and orientation disposed to deliver an induced 
electric field substantially perpendicular to said cerebral cortical region of 
interest; 



positioning said subject to immobilize said cerebral cortical region of interest; 



5 



directing said stimulator to said subject so that said stimulator is focused to said 
cerebral cortical region of interest, said directing based on said determined 
stimulator position and orientation; and 



delivering said TMS to said subject, said delivering comprising generating an 

electric field substantially perpendicular to said cerebral cortical region of 



10 



interest. 



19. The method of claim 1 8, further comprising performing at least one spatial 
calibration of a delivery device to a data set of said subject. 

20. The method of claim 1 9, wherein said directing is further based on said at 
least one spatial calibration. 



comprises a digitizer to a subject image calibration. 

22. The method of claim 1 8, wherein said defining comprises using a 
functional image to obtain said region of interest. 

23. The method of claim 1 8, wherein said determining comprises calculating 
20 said stimulator position and orientation based on a scalar product of a model of said 

induced electric field times a vector model of said region of interest. 

24. The method of claim 22, wherein said determining comprises 
superimposing a model of said stimulator on a model of said functional image and an 
anatomical image. 

25 25. The method of claim 24, wherein said directing further comprises 

modeling said immobilized subject to obtain a scalp surface model and registering said 
scalp surface model to said model of said functional image and said anatomical image. 



15 



21 . The method of claim 19, wherein said at least one spatial calibration 
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26. The method of claim 18, wherein said directing is performed by a robotic 
arm to which said stimulator is adapted. 

27. A method of planning optimal delivery of TMS, comprising: 

imaging a subject to obtain an anatomical image, said anatomical image including 
5 a representation of a region of interest; 

obtaining a functional image of said subject in which said region of interest is 

selectively activated; 
identifying a target point from said region of interest via said functional image; 
modeling said anatomical image to obtain an anatomical model; 
10 co-registering said target point with said functional image and said anatomical 

model; and 

determining a position and orientation of a TMS stimulator so that maximum 
biological efficacy of said target point is obtained. 

28. The method of claim 27, further comprising storing said position and 
15 orientation. 

29. The method of claim 27, further comprising providing said position and 
orientation to a computer, said computer for control of a robot to which TMS stimulator 
is adapted. 

30. The method of claim 29, further comprising moving said robot into a 
20 location such that said TMS stimulator is placed in said position and orientation. 

3 1 . The method of claim 30, further comprising delivering said TMS. 

32. The method of claim 3 1 , wherein said TMS delivery comprises 
preoperative mapping. 

33. The method of claim 31, wherein said TMS is delivered during imaging. 

25 34. The method of claim 3 1 , wherein said TMS is delivered for treatment of a 

specified condition. 
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35. The method of claim 34, wherein said specified condition comprises long 
term depression. 

36. The method of claim 31 , wherein said TMS delivery comprises a plurality 
of individual treatment sessions over a plurality of days. 

5 37. A magnetic stimulation device, comprising: 

a first wing having a first plurality of windings, said first plurality of windings 
having first portions and second portions, said first portions having a 
smaller spacing therebetween than said second portions; 
a second wing having a second plurality of windings, said second plurality of 
10 windings having third portions and fourth portions, said third portions 

having a smaller spacing therebetween than said fourth portions; and 
wherein said first wing is adjacent said second wing such that said first portions 
are adjacent said third portions. 

38. The magnetic stimulation device of claim 37, wherein said magnetic 
1 5 stimulation device comprises a minimum inductance. 

39. The magnetic stimulation device of claim 38, wherein said magnetic 
stimulation device comprises a minimum power dissipation. 

40. A magnetic stimulation device, comprising: 

a two-winged coil having a plurality of wires, said plurality of wires being 
20 concentrated at a center portion of said two-winged coil and extending 

from said center portion in a plurality of arcs. 

4 1 . The magnetic stimulation device of claim 40, wherein said plurality of arcs 
extend from said center portion at an increased spacing therebetween. 

42. The magnetic stimulation device of claim 40, wherein said plurality of 
25 wires comprises a first portion and a second portion, said first portion being adapted to 

cany current in a direction opposite to said second portion. 
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43. A magnetic stimulation device, comprising: 

a coil comprising a first plurality of turns having a positive polarity and a second 
plurality of turns having a negative polarity. 

44. A magnetic stimulation device, comprising: 

a first layer of wires adapted to allow a first current to travel therethrough; and 
a second layer of wires adapted to allow a second current to travel therethrough, 
said second current being opposite to said first current. 

45. A magnetic stimulation device, comprising: 

a coil having a first wing, a second wing, a third wing, and a fourth wing, each of 
said wings extending peripherally from a center portion of said coil; and 

wherein each of said first wing, said second wing, said third wing, and said fourth 
wing comprises a proximal portion adjacent said center portion and a 
distal portion, said proximal portion having a first plurality of windings 
having a first spacing therebetween, said distal portion having a second 
plurality of windings having a second spacing therebetween, said second 
spacing being greater than said first spacing. 

46. The magnetic stimulation device of claim 45, wherein each of said first 
and second plurality of windings are generally concentric. 

47. A method of designing a magnetic stimulation device, comprising: 
defining a desired electric field at a preselected point in a volume, said preselected 

point corresponding to a first portion of said magnetic stimulation device; 
and 

arranging a plurality of windings to obtain said desired electric field at said 
preselected point. 

48. A magnetic stimulation device, comprising: 

a first central coil and a second central coil, each of said first central coil and said 
second central coil extending outwardly from a proximal portion to a distal 
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portion, said proximal portion of said first central coil adjacent said 
proximal portion of said second central coil; and 
a plurality of outer coils adapted around said first central coil and said second 

central coil, each of said plurality of outer coils extending outwardly from 
a proximal portion to a distal portion, said proximal portion of each of said 
plurality of outer coils adjacent said proximal portion of said first and 
second central coils. 

49. The magnetic stimulation device of claim 48, wherein said plurality of 
outer coils comprise four coils. 

50. The magnetic stimulation device of claim 48, wherein said first central coil 
and said second central coil comprise positive turns and negative turns. 
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